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Abstract 

Recent experiments indicate a departure from the exact Tri-Bimaximal mixing by measure ring definitive non-zero 
value of #13. Within the framework of type I seesaw mechanism, we reconstruct the triangular Dirac neutrino 
mass matrix from the fi — T symmetric mass matrix. The deviation from fi — t symmetry is then parametrized by 
adding dimensionless parameters yi in the triangular mass matrix. In this parametrization of the neutrino mass 
matrix, the non-zero value #13 is controlled by Ay = 2/4 — yg . We also calculate the resulting leptegenesis and 
show that the triangular texture can generate the observed baryon asymmetry in the universe via leptogenesis 
scenario. 
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1 Introduction 

The discovery of the neutrino masses and the lepton mixing have played a pivotal role in probing physics beyond 
the Standard Model (SM). Recent experiments MINOS [U|2], T2K [3], Double CHOOZ @], Daya Bay [5] and 
RENO [5] have reported definitive non-zero #13 result. The latest global analysis of neutrino oscillation data 
yields the following best-fit values with la, 2a and 3a errors for the oscillation parameters [7]. 
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Table 1: Global oscillation analysis with best fit for Am^, Am^, sin 2 #12, sin 2 6*23, sin 2 6*13 and S the upper and/or 
lower corresponds to normal and/or inverted neutrino mass hierarchy. 
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Seesaw mechanism [5] gives a natural explanation of the smallness of the masses for light neutrinos and mixing 
by connecting the tiny neutrino masses to a very heavy right-handed neutrinos masses. 

According to the seesaw mechanism, lepton number is broken at high energies due to right-handed neutrino Ma- 
jorana masses, resulting in small left-handed neutrino Majorana masses suppressed by the heavy mass scale. The 
seesaw mechanism also provides an attractive mechanism for generating the baryon asymmetry of the universe 
via leptogenesis HH] ■ 

The seesaw Lagrangian can be written as : 

C = IrMiIl + NrMdvl + ^NrMrNr + h.c. (1) 

where vl = (z/ e , z/^, z/ r )'^ 1 , II = (e, /j,t) t and Nr = (iVi, 7V 2 , -/V 3 ) denotes the left-handed (light) neutrinos, the 
left-handed charged leptons and the right-handed (heavy) Majorana neutrinos, respectively. We assume that 
the mass matrices of both the heavy Mr = diag(M\, M 2 , M 3 ) and the charged lepton Mi are diagonal and real 
matrices. 

After integrating out the heavy right-handed neutrinos, the symmetric Majorana mass matrix for light neutrinos 
is : 

M v = -MdM^MZ. (2) 
The mixing matrix V — UPuaj that diagonalizes M u is : 

M v = VMf ag V T (3) 

where Af*° 9 = diag(nii, 7712,7713) and 

C12C13 S12C13 si 3 e~ l<5 

U = -S12C23 - c 12 s 23 s 13 e lS c 12 c 23 - s 12 s 23 s 13 e lS s 23 c 13 (4) 
S12S23 - c 12 c 23 s 13 e lS -C12S23 - si 2 c 2 3Si 3 e lS c 23 ci 3 

where = cos 6ij , = sin 9^ , 8 is the CP-violating phase and Puaj is diagonal matrix consisting of non-trivial 
Majorana phases. 

Now any nonsingular matrix can be decomposed into the product of a nonsingular, upper or lower triangular 
matrix Tljj and a unitary matrix W. In particular, the Dirac matrix Md can be written as : 

M D = WT LtU . (5) 

It has been noticed that the relevant quantities for leptogenesis are the CP-violating phases in M^M d = T^ v Tljj 
which means that the unitary matrix W cancels out in the unflavored leptogenesis but not in the seesaw formula. 
Our ansatz consists in taking the unitary matrix W to be identity implying that the Dirac matrix is cither lower 
or upper triangular matrix Mu = Tljj. 

Triangular textures have been considered as the simplest form to study since all the unphysical features can be 
eliminated from the start [TT] and [H] . Leptogenesis with triangular ansatz has been studied by [131 Q3] . 

The paper is organized as follows. In the next section, we consider a lower or upper triangular matrix as Dirac 
mass matrix. We show how to reconstruct the triangular Dirac matrix from the neutrino symmetric mass matrix 
for the Tri-Bimaximal (TBM) mixing. In section 3, we discuss the breaking of the triangular mass by adding 
dimensionless parameters. We next analyze the mass texture by fitting it to the observed neutrino oscillation 
data. We show that the triangular mass matrix can naturally accommodate the observed neutrino oscillation 
parameters. In section 4, we calculate the relevant quantities for leptogenesis and estimate the baryon asymmetry 
that is consistent with the observed value. 
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2 Triangular Dirac Matrix for Tri-Bimaximal Mixing 

In the following, we reconstruct the Dirac matrix Mp = T® u analytically from the neutrino symmetric mass 
matrix elements rriij ( i.e. light neutrinos masses, heavy Majorana masses and the Majorana phases ) for the 
Tri-Bimaximal mixing pattern #23 = j , #12 = sin -1 (l/y/3) and #13 = [15] ■ The general symmetric matrix that 
is diagonalized by the Tri-Bimaximal mixing matrix V° can be written as : 



M?, 



-T r ° rrM^T, 




(6) 



with the eigenvalues 



uin = a 



2b 



2c — a — b. 



(7) 



By solving the above equation, we get the following Dirac matrix for upper or lower triangular matrix : 
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(8) 



where d! = det(M v ) is the determinant of the neutrino mass matrix. 



3 Deviations from Tri-Bimaximal 

Now we consider deviation from TBM by breaking the triangular neutrino matrix with dimcnsionless parameters 
Vi as : 
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(9) 



with a' = a e z<Pa ,b' = b e l<Pb and d = c 
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It is interesting to notice that by taking j/3 = yi = 2/6 and 2/5 = 2/4 for the lower triangular and 2/3 = 2/2 = 2/1 
and 2/5 = 2/4 f° r upper triangular, the neutrino matrix M v will be : 

(a' y\ b' 2/12/4 b' 2/12/6 \ 

6' 2/12/4 d y\ {a' +b'- d) 2/42/6 • (10) 

V 2/12/6 («' + &' - c') ?/42/6 c' y\ J 

For 2/4 = 2/6 = 1, M„ has \x — r symmetry leading to #13 = and #23 = f • Moreover, in the limit 2/1,2/4,2/6 — > 1, 
the above matrix gives the TBM angles and mass eigenvalues. 

In particular for 2/1=2/2=2/3=2/4 = 2/5 = 1, the lower or upper triangular mass matrix gives the neutrino matrix 
M u : 

( a' b' b'y 6 \ 

M v = [if d (a' + b'-d)y e (11) 

V b' y 6 (a' + b'- d) 2/6 d 2/1 / 

which is the strong scaling Ansatz studied by [16]. 

To see how neutrino mass matrix given by Tljj can lead to deviations of neutrino angles from their TBM 
values, we consider the hermitian matrix H v = M v Ml, 

( Vi A 2/12/4-B 2/12/eC \ 
H v = ym B* y\D y 4 y 6 E (12) 
V 2/12/eC* 2/42/e^* ) 

where 

A = y 2 a 2 + (yl + y 2 6 )b 2 

B = yia'b'* +ylb'd* + ytb\a'* + b'* -d*) 

C = y\o!b'* + ytb'{a'* + b'* -d*) + y 2 b'd* 

D = y 2 b 2 + y 2 c 2 +y 2 \a' + b'-c'\ 2 (13) 
E = y 2 b 2 + y 2 c\a'*+b'*-d*) + y 2 {a' + b'-d)d* 

F = y 2 b 2 + y 2 \a' + b>-d\ 2 +y 2 c 2 . (14) 
As a consequence, the hermitian matrix B v is diagonalized by the unitary matrix V, 

H v = vuf as Mf as tyt (15) 
where the deviation of the Oij from their TBM values is parametrized by three quantities = 1, 2, 3) : 

6*23 = ±^ + £ 2 3 #13 = £13 #12 = sin -1 (~^J + e 12 ■ (16) 

The angle #23 of the atmospheric mixing is expressed in terms of the parameters 2/1 as : 

tan023 - Mcj 

2/4 a (2/1 - 2/1 ) sin fab ~ c(yj - y%) sin ip cb 
2/6 a(y 2 ~ yj) sin (p ab + c(yj - y%) sin <p cb 

where ip ab = ¥>a - <Pb and ip cb — ip c - ip b . 



(17) 
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The Dirac CP phase S is given by : 

2/42/6 Im(E) 



tan S 



C23S23(yl D ~ vl F ) + 2/42/6 (c| 3 - s 2 ^)Re{E) 
2cy i y 6 (yl - y 2 ) [a sin ip ab - 6sin^ c6 ] 



where 



The atmospheric angle #23 is 



2/6 

and the deviation 623 of the atmospheric mixing will be : 



We see that the deviation £23 vanishes for j/4 = y§ as it should be. 
The Dirac CP phase 8 becomes : 



c{vl +2/g)siny? e 
ayj + c(yj + y%) cos <p c 



A useful measure of CP violation is given by [17) : 



J cp = - sin 2^12 sin 2^23 sin 6*i3 cos 6>i3 sin (5 
8 

Im(H v \ 2 H V 23-Hy 31) 



tai12 ' 13 - a 2/ 1 (2/? + 2/4 2 + 2/6 2 )V Vi^yl • ( } 



The solar mixing angle 6*12 is governed by : 



22/12/6 VX 

a{yi + vl + vD^vlvl - 2/12/4 - vlvl) 



(18) 



A = (yl-y 2 )[b 2 y 2 + c 2 (y 2 +y 2 )]sm(2e 23 ) 

+ 2y 4 y 6 [b 2 yj + c(yl+yl)(~c + asmtp ab + bcosip cb )]cos(2e 23 ) ■ (19) 

In this analysis, we consider only the case of a' ~ b' . In the case, we have Im(C) ~ —Im{B). 



tan 023 - -— (20) 



, 2/4-2/6 , 01 V 

tan e 2 3 — ■ (21) 

2/4 + 2/6 



(23) 



Aniji Amjj Am§ 2 
Here we have : 

Im{H vl2 H u2Z H vZl ) - 2a 3 cylylyl{yl - yl){ylyl + y\yl + Aylyl){a 2 +c 2 -2ac cos ip ca ) smip ca (24) 
which shows that for j/4 = y§, the rephasing invariant Jcp vanishes. 
The reactor mixing angle #13 is : 



tan 26»i2 ^ „/„,2 i ..2 i 3T7T^72 IX^ ( 26 ) 
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where 

X = c(yl-y%) 2 (yl+yl)[c(yi + yl)+2ayl cos ip ca ] 

+ a 2 [4yl(yj + y 6 2 ) 2 (2y 2 + y\ + y 6 2 ) + y*(5yi + 2y 4 2 y 6 2 + y 6 4 )] (27) 

which for y\ = j/4 = ye, we get tan26'i2 = 2\/2 as expected. 

Finally, the mass squared differences are given by : 

Am 2 21 ~ a 2 (y 2 + 2y 2 )[y 2 + 2y 2 + 4y 6 (y 4 -y 6 )] +0((y 4 -y 6 ) 2 ) 

Am 2 31 ~ 4^( c 2 +a 2 -2accos^ ca )[y 6 + 2( 2 /4-2/ 6 )] + 0(( 2 /4-y6) 2 ) • (28) 

In the numerical analysis, we choose the parameters yt, a, c and tp ca as inputs. With the help of the 3cr allowed 
experimental neutrino results, we determine the allowed ranges of these inputs. 

In Fig.l, we have plotted €23 with respect to ye and Ay = 7/4 — y$ and obtain restriction on the parameter space 
of ye and Ay. One can see that ye and Ay are constrained in the ranges 0.97 < ye < 1.01 and —0.2 < Ay < 0.2. 




Figure 1: Deviation of the atmospheric angle €23 versus ye and Ay. 

In the limit y\, 7/4, ye — > 1 and a' ~ 6' <C c', the mass-squared differences Am 2 ^ and Ato^ are : 

Am^! ~ 9a 2 

Am 2 ! ~ 4c 2 . (29) 
Now, using the best fit of the masses-squared differences, we get : 

a ~ 2.91 x 10~ 3 eV 

c ~ 2.51 x 10~ 2 eV. (30) 

In Fig. 2 and Fig. 3, we display the exact expressions for the mass-squared differences versus Ay and the phase 
ip C a for yi = 0.9 and ye = 0.99. It can be seen that for these values Ato 2 i an( l Ato 2 ^ are within 3<r ranges. Since 
Am§ x depends on the cos<y9 CQ , Fig. 3 shows two symmetrically regions by varying ip ca between and 2~k. 
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r 



For ?/i = 0.9 and ye = 0.99, we show in Fig. 4 and Fig. 5 the variation of the reactor #13 and atmospheric 
ei2 angles versus Ay and (p ca . These plots indicate that $13 and €12 are within 3cr ranges. Moreover, since #13 
depends on the Ay, Fig. 4 shows two symmetrically separated regions within these ranges. 
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4 Leptogenesis 

The CP violating asymmetry is due to the interference of tree level contribution with the one-loop corrections for 
the decay of the ith heavy neutrino HH] : 

e ' 1 r(jvi->#<o + r(ivi->0tj) ■ 1 > 

The contribution of the heavier neutrinos is washed out and only the asymmetry generated by the lightest neutrino 
survives : 

Here v ~ 174 GeV is the electroweak symmetry breaking scale and the function / stems from the contributions 
of both self-energy and vertex diagrams : 

fix) = ( T ^ + l-(l-x)Ml + W) ■ (33) 

For hierarchical heavy Majorana neutrinos i.e., Mi -C M2,3, one has : 

/W - -jjg ■ (34) 

The CP asymmetry ei gives rise to a lepton asymmetry Yl in the universe : 

n-L — TIt K 

Yl = — = — ei 35 

where s denotes the entropy density, g± ~ 106.75 is the effective number of relativistic degrees of freedom 
contributing to the entropy of the early universe and k is the dilution factor that is obtained by solving the 
Boltzmann equations. To a good approximation, n is given by |18j : 

1 jqgfr* 10<^< 10 6 (36) 

where the parameter K, which is defined as the ratio of the thermal average of the N\ decay rate and the Hubble 
parameter at the temperature T = Mi, is : 

K = MP 2 V&Vnhi . (37) 

1.66 x Sirv 2 ^ Mi v ' 

Here M P ~ 1.22 x 10 19 GeV is the Planck mass. 

The produced lepton asymmetry Yl is converted into a net baryon asymmetry Yb through the non-perturbative 
sphalcron processes. The baryon asymmetry is obtained : 

£ 8N f + 4Nh 

where in the standard model Nf = 3 and Njj = 1 are the number of fermion families and complex Higgs doublets 
respectively. Typically, one gets : 

00 

Yb^- T9 Yl ■ (39) 



9 



For the lower triangular Dirac neutrino mass matrix with a ~ b -C c, the CP asymmetry is 

3cMi(i£-yg) 2 



ei 



The value of -R' for the lower triangular texture is expressed as : 



sin ip c 



K 
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/ 2 , 2 i 2\ 
(2/1 +2/4+%) 



(40) 



(41) 



which is independent of M\ and is approximately equal to 11.3 
The baryon asymmetry generated by leptogenesis is : 



3.9 x 10- 6 (yf - yl) 2 M : 



sin ip c 



(42) 



(y? + y 2 + % 2 ) 2 (1-3 + in(y 2 + 2/I + % 2 )) - 6 

We see that for an exact y, — r symmetry, the baryon asymmetry 77 s vanishes |19j . 

We use r\B = (6.19 ± 0.15) x 10~ 10 as the upper and lower bound of the baryon-photon ratio from the WMPA 
observation [50]. Our result is shown in Fig. 6 with Mi = O(10 13 ) GeV, where the two horizontal lines represents 
the WMPA allowed bounds. We can see that the triangular Dirac mass matrix provides a baryon asymmetry 
consistent with the current observations. 




Figure 6: Baryon asymmetry r\B versus ip ca for Mi = {1.5, 2, 2.5, 3} x 10 GeV and Ay — 0.08 



5 Conclusion 

In this paper, we have reconstructed lower and triangular textures of the neutrino Dirac mass matrix that can be 
obtained via type I seesaw mechanism from the neutrino symmetric mass matrix for the TBM mixing pattern. 
The recent experiments have reported definitive non-zero value for the reactor mixing angle. Therefore, we have 
considered deviation of the triangular ansatz from TBM by adding dimensionless parameters. We have studied 
the phenomenological implications of these triangular textures and have obtained interesting results for the mixing 
angles and mass-squared differences. 

Furthermore, we have shown that the breaking of the triangular mass matrix from the exact TBM gives a non-zero 
baryon asymmetry via leptogenesis consistent with the current observations. 
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